Observation and characterization of laser-driven 
Phase Space Electron Holes 



o 
o 

(N 

< 

m 



<d : 

C/3 ! 

C ■ 

• i-H , 

O 1 

• i-H ' 

CZ5 ■ 
>Y 

^ : 

Ph. 



> 
(N 

in 

od 
O 

o 



G. Sarri 1 , M. E. Dieckmann 2 , C. R. D. Brown 3 , C. A. Cecchetti 1 , D. J. Hoarty 3 , 
S. F. James 3 , R. Jung 4 , I. Kourakis 1 , H. Schamel 5 , O. Willi 4 and M. Borghesi 1 

1 School of Mathematics and Physics, The Queen's University of Belfast, Belfast, BT7 INN, UK 
2 ITN, Linkoping University, 60174 Norrkoping, Sweden 
3 AWE, Aldermaston, Reading, Berkshire RG7 4PR, UK 
4 Institute for Laser and Plasma Physics, Heinrich- Heine- University, Diisseldorf, Germany 
5 Physikalisches Institut, Universitdt Bayreuth, D-95440 Bayreuth, Germany 

The direct observation and full characterization of a Phase Space Electron Hole (EH) generated 
by laser-matter interaction is presented. This structure has been detected via proton radiography 
during the interaction between an intense laser pulse (r ~lns temporally flat-top, 1~ 10 14 W/cm 2 ) 
and a gold 26/xm thick hohlraum. This technique has allowed us the simultaneous detection of 
propagation velocity, potential and electron density spatial profile across the EH with fine spatial 
and temporal resolution providing an unprecedentedly detailed experimental characterization. 



Phase-Space Electron Holes (EH) [l[ are electrostatic 
excitations in collisionless plasmas characterized by a 
positive potential hump in which a population of elec- 
trons is trapped. In addition to their relevance to many 
fundamental plasma processes such as two-stream insta- 
bilities [2| and saturation in Landau damping [3j , EH play 
a key role in a wide range of space plasma scenarios (e.g. 
microscopic dissipation during magnetic reconnection in 
the Earth's magnetosphere [J] or the generation of cos- 
mic ray electrons in supernovae [5|) and are commonly 
detected in near Earth plasmas [y, 0|- The detrimen- 
tal effect they have on the focusing properties of parti- 
cle accelerators and storage rings has also been recently 
highlighted [§']. This omnipresence of such structures 
in collisionless plasmas requires therefore, beside a deep 
theoretical study [![, an equally profound experimental 
characterization. Previous experimental work detected 
this type of structures in magnetized collisionless plas- 
mas during either strong discharge in gas tubes or 
magnetic reconnection in toroidal plasma current sheets 
[Io| . In all these cases the existence of EH was deduced 
from positive spikes in high-bandwidth Langmuir probes. 
Remarkably advancing these previous detections, here we 
report the first direct observation of an EH in a laser- 
matter interaction experiment, and suggest a new way to 
generate and study them in a controllable manner. EH 
diagnosis using proton radiography [Til . Il2| allows in fact 
the simultaneous measurement of propagation velocity, 
potential and electron density spatial profile and tempo- 
ral evolution leading to a comprehensive characterization 
with fine temporal and spatial resolution, of the order of 
few picoseconds and few microns respectively 

The present experiment was performed at the He- 
len laser system in AWE [l|. It involved the illumi- 
nation of the inner surface of a hohlraum target by an 
intense and relatively long (r wins temporally flat-top, 
I« 10 14 W/cm 2 , A=0.527nm) laser pulse, an interaction 
which produces conditions relevant to Indirect Drive In- 
ertial Confinement Fusion (ICF) experiments [l4j| . The 
hohlraum target consisted of an open-ended Au cylinder 
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Figure 1: a. Experimental setup: The experiment, car- 
ried out in the Helen facility [L§], used two different laser 
beams. The EAST beam (I w 10 14 W/cm 2 , r w Ins flat- 
top, A ~ 0.527/im) irradiated the inner surface of a 1.5mm 
diameter, 26/im-thick open-ended gold hohlraum. The in- 
teraction was probed by the proton beam generated in the 
interaction between a gold foil and the CPA beam 

(I > 10 19 W/cm 2 , r » 700fs). b.Proton Radiography of 
unirradiated hohlraum: the Laser Entrance Hole (LEH) 
is visible at the upper-left side of the hohlraum. 



with diameter of 1.5 mm, length of 1 mm, and wall thick- 
ness of 26 /zm. The interaction was probed via proton ra- 
diography [111, E3] and a brief sketch of the experimental 
setup is visible in Figure[T] The interaction beam (EAST 
beam in Figure La) was focused through a 350/im diame- 
ter entrance hole onto the inner surface of the hohlraum. 
A second short and intense pulse (r w 700fs, I> 10 19 
W/cm 2 , CPA beam in Figure La) was focused onto a 
20/im gold foil in order to create, via Target Normal 
Sheath Acceleration (TNSA), [l5| a wide spectrum pro- 
ton beam, with a changeable delay with respect to the 
EAST beam. The proton beam, after having probed the 
plasma, was recorded on a stack of RadioChromic Films 
(RCF), [l6j]. Such a probing scheme enables the moni- 
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Figure 2: a. Experimental data: RCF image associated to a time ~ 160ps after the beginning of interaction. The hohlraum 
is open-ended and has a diameter of 1500/im. The dashed square in the top-right corner highlights the feature discussed in this 
Letter. The white ellipses outline zones of proton density depletion due to the electrostatic charging of the hohlraum walls. The 
white ruler represents 500/im. b. Zoom of the structure: Zoom of the region outlined by the dashed square in a). The white 
ruler represents 100/xm. c. Proton signal: modulation in proton dose on the RCF layer in correspondence of the structure. 



toring of the transverse electric field distribution inside 
the plasma by measuring the deflection of a proton beam 
as it passes through it. The technique exploits the fact 
that, as a consequence of the high degree of laminarity 
of the beam, the proton source, while being physically 
extended, is practically equivalent to a nearly point-like 
virtual source located in proximity to the target [13 • A 
point projection of the probed region is thus obtained 
with a geometrical magnification of M ~ 10. Under the 
assumption of small deviations (i.e. the proton trajecto- 
ries do not cross) , the transverse electric field distribution 
can be derived directly from the relative modulation of 
the proton density deposited on a given RCF layer pd| : 
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where < Ey > is the transverse electric field compo- 
nent averaged along the longitudinal dimension, e p is the 
probe proton energy, L is the distance between the inter- 
action area and the detector, b is the longitudinal length 
of the non-zero electric field region and dn p /n p is the rel- 
ative modulation of the proton density at the detector 
plane. 

Data exemplifying the features observed by proton 
imaging are displayed in Figure [51 a and [21 b. As a general 
rule, the electric fields are directed from the regions of a 
lighter grey color compared to the background (reduced 
probe proton flux) towards the regions of darker grey 
color (increased flux). A region of pronounced modula- 
tion in the probe proton density (Figure HJc), evidence of 
a modulated electric field distribution, is observed ~300- 
400/^m from the rear surface of the irradiated target sur- 
face. A region of turbulent plasma, resulting in a chaotic 
deflection pattern, is visible at the interaction point even 
though the aim of the Letter is to study the solitary 
structure at the rear surface of interaction. Thanks 
to the broad spectrum of such a proton beam, differ- 
ent layers within the RCF stack provide snapshots of 
the interaction at different times even in a single shot 
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Figure 3: EH propagation as a function of time relative 
to the beginning of the interaction. The temporal error bars 
are due to the finite transit time of the protons across the 
structure. Experimental data indicate a constant velocity of 
(1.6 ± 0.6) ■ 10 6 m/s~ (1.6 ± O.Q)v th . 



configuration, 11] . This density modulation (shown in 



Figure Oc) shows to propagate with a constant velocity 
of v ~ (1.6 ±0.6) • 10 6 m/s (Figure[3]) while maintaining a 
substantially time-independent profile in the co-moving 
reference frame. The electric field distribution across the 
structure, E(x), has been extracted (Figure [Ha) using 
Equation[T] The corresponding potential profile has been 
calculated by spatial integration of E(x) (Figure [Hb). In 
this calculation we have assumed a quasi-planar structure 
with longitudinal dimension b of the order of the trans- 
verse dimension, i.e. « 600/im; these symmetry con- 
siderations concur with published numerical results [l8l |. 
The potential exhibits a localized bell-shaped structure 
80 — 90/im wide with a maximum value of « 30V. 

In order to estimate the plasma parameters in the re- 
gion of observation, the interaction between the laser 
and a 26/im thick gold foil has been simulated using a 
ID hydrodynamic Lagrangian code (Hyades) including 
radiation transport and ionisation [191 ]. Simulations in- 
dicate that the more energetic X-rays generated during 
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Figure 4: a) Electric field distribution: experimental 
electric field distribution across the soliton structure (blue 
empty circles) compared with the electric field profile (red 
line) obtained, by using Poisson's equation, from the theo- 
retical potential in Figure b). b) Potential spatial pro- 
file: experimental potential profile (blue empty circles) com- 
pared with the theoretical potential profile (red line) obtained 
by adding the contribution of two partially overlapping EHs 
(green dashed curves), c) Electron density depletion: ex- 
perimental electron density distribution inside the EHs; two 
different regions of electron depletion are clearly visible sug- 
gesting the presence of two different EH structures. 



the interaction propagate through the gold foil and ion- 
ize the ambient gas at the rear surface (pressure ~ 1CP 3 
mbar) creating a steady plasma in thermal equilibrium; 
the electron temperature and density are predicted to 
be n e ~ 2.5 • 10 12 cm -3 and T e ~ 2eV respectively, 
implying a Debye length of Xd — 7/zm, an electron 
plasma frequency of u> pe ~ 10 11 s _1 (ion plasma frequency 
Lo V i ~ 5-10 8 s _1 ), an electron thermal velocity of vth — 10 6 
m/s and an ion acoustic velocity of c s ~ 3 • 10 3 m/s. 
The velocity and the width of the structure are then 
~ l.dvth ( or — 530c s ) and ~ 10 — 12A_d respectively, 



while the normalized maximum value of the potential is 
(f) = eV/KsTe ~15. According to these parameters, the 
ratio of the electron mean free path to the electron Debye 
length is approximately 3 • 10 3 . The plasma is then colli- 
sionless and it can thus support propagating EHs .IS]. 
Since the plasma is probed at ss 100 — 200ps after the 
beginning of interaction and the ion plasma frequency is 
ujpi ~ 5 • 10 8 s _1 (w^f ~ 2ns), it is reasonable to neglect 
motion of the ions and to consider them as a still, neutral- 
izing background in the plasma. Under this assumption, 
the electron density across the structure can be extracted 
from the data; the charge density, obtainable by Poisson's 
equation, is equal to p(x) — e ■ (r^(ac) — n e (x)) where 
n e [x) (rii(x)) is the electron (ion) density distribution. 
Simulations indicate an ionization state Z — 1 therefore 
the ion density within the structure can be expressed as: 
rii{x) = n i0 = n e0 leading to n e (x) = n i0 - p(x)/e, (Fig- 
ure IHc). The electron density exhibits two pronounced 
depleted regions suggesting the simultaneous presence of 
two partially overlapping EHs within the structure. The 
electron depletion is of the order of ss5-10%, well within 
the range of density depletions seen in published simula- 
tions ranging from 1-2% [lj] up to 15-20% [loj]. 

For EHs in unmagnetized plasmas, the potential pro- 
file can be described in a first approximation by its small 
amplitude expression: 

<j>(x) = ^maxsech 4 ( — J , (2) 

where j e is a numerical factor ranging from to 1 
depending on the electron velocity distribution in the 
plasma |21| . j e =1 corresponds to a pure Maxwellian 
distribution while a smaller value of 7 e implies a big- 
ger deviation from the pure Maxwellian behavior or, 
equivalently, a bigger value of k in the K-distribution 
[22j |. The experimental potential shape (Figure HJb) 
can then be interpreted as the sum of two different 
bell potentials with 7 e w 0.7 corresponding to k ~ 4; a 
strong deviation from a pure Maxwellian distribution 
can reasonably be due to the hot-electron population 
created during the laser interaction with the hohlraum 
walls. The comparison between the experimental and 
theoretical profile is shown in Figure [4jb where the two 
single potential profiles are also plotted. Differentiating 
this theoretical potential profile leads to a theoretical 
electric field distribution that reproduces the behavior of 
the experimental one (Figure [4] a). This interpretation 
is corroborated by simulations reported in literature [l8| 
showing indeed the creation of 'trains' of EHs when an 
electrostatic perturbation is applied at the edge of the 
simulation box. The maximum value of the velocity 
distribution of the electrons trapped by the EH can be 
roughly estimated by equating the maximum value of 
the potential energy with the electron kinetic energy: 
v te = (2eV max /m e ) 1 / 2 ~ 6.7 • 10 6 m/s~ 6.7v th . The 
velocity oscillations within the electrostatic potential ex- 
ceed the observed propagation speed (1.6 ± 0.6) • 10 6 m/s 
of the structure by a factor 4, implying that this 
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potential must trap the electrons of the ambient plasma, 
which is at rest. The velocity oscillations of the trapped 
electrons exceed the simulated thermal speed by a factor 
6—7, which is in accordance with theory [l[ and also 
not unusual for an EH in simulations 0]. 

The perturbation driving the EH generation is thought 
to be the sudden charging of the hohlraum walls due 
to the residual positive charge in the target left behind 
by the accelerated electrons that are energetic enough 
to escape [H, |23 |. This potential can be estimated 
from the deflection of the protons passing close to the 
walls (highlighted in Figure [5] a). Contiguous regions of 
proton depletion and accumulation are in fact present 
in proximity of the walls consistent with an electrostatic 
potential at the wall surface of w 800V = AQQK B T e /e. A 
significant difference in amplitude between the EH and 
the driving potential is reported both in experiments 
and simulations [l8[ and it thus appears to be a 
necessary condition for the excitation of such structures. 

The experimental propagation velocity of (1.6 ± 0.6)vth 
slightly exceeds the allowed velocity range 
veh < 1.307v t /j set by the analytical theory devel- 
oped by Schamel et al. [1(. However, it has to be noted 
that this range is valid only for a pure Maxwellian 
distribution and for EHs that have reached stationarity. 
The significant deviation from a Maxwellian distribution 
that the data suggest, together with the very early 
times at which the structure is observed (t ■ oj pe « 16), 
could be a possible explanation for this partial disagree- 
ment. Indeed, simulations by Califano et al. [l8| show 
propagation velocities ranging from a fraction of Vth up 
to 2v t h- The extremely high ion Mach number of the 
structure ( M=530) definitely excludes the possibility of 
an ion wave driven phenomenon such as a collisionless 
shock [H]. 

Published analytical models [l[ and numerical simu- 
lations [TH, [2(| set the conditions to identify a solitary 
structure in a plasma as an EH to be: a positive hump- 
like potential, a depletion in the electron density at the 
potential peak and finally a propagation velocity of the 
order of the thermal electron velocity. The results shown 



in Figures [3] and [4] satisfy all the requirements expressed 
above. Even the spatial width of « 10 — I2\rj concurs 
with the hypothesis of an EH detection; simulations 
[HI and space plasma measurements Q show that these 
structures have a typical width ranging from a few up 
to tens of Debye lengths. 

In summary, we have reported the first laboratory ob- 
servation of EHs in a laser-plasma experiment. Thanks 
to its fine temporal and spatial resolution, of the order 
of few picoseconds and few microns respectively, proton 
radiography allows for a comprehensive characterization 
of the structure. Data analysis indicates the generation 
of two partially overlapping EHs moving with a velocity 
close to the thermal electron velocity. Electric field, 
potential profile and electron density distributions have 
been measured, together with the propagation velocity 
and width of the structure, giving results in good agree- 
ment with published analytical theory and simulations. 
Temporal evolution of the structure has also been 
accurately measured. The experimental setup shown 
permits to create and comprehensively characterize 
EHs, or other soliton-like structures, that have many 
implications in a wide range of plasma scenarios. Such 
structures are one of the most widespread nonlinear 
phenomena in collisionless plasma ranging from Earth's 
vicinity and supernova scenarios to particle accelerator 
storage rings. By controlling the characteristics of the 
ambient plasma or the laser parameters, it may be 
possible to study this type of phenomena in a wide 
range of different regimes opening up the possibility of 
acquiring unprecedented experimental knowledge in this 
field. 

Funding for this research has been provided by 
the AWE Academic Access Scheme, EPSRC grants 
EP/E035728/1 and EP/C003586/1, by DFG TR 18, 
GK 1203 and FOR 1048 VR. The authors acknowledge 
the support of the HELEN Facility laser and target 
preparation personnel and K.Quinn for the fruitful help 
given. 



[1] H. Schamel Phys. Rep. 140,161 (1986). 

[2] M. E. Dieckmann et al. New J. Phys. 8,225 (2006). 

[3] J. R. Danielson et al. Phys. Rev. Lett. 92,245003 (2004). 

[4] J. F. Drake et al. Science 299,873 (2003). 

[5] K. Koyama et al. Nature 378,255 (1995). 

[6] F. S. Mozer et al. Phys. Rev. Lett. 79,7 (1997). 

[7] M. Hoshino Science 299,834 (2003). 

[8] M. Blaskiewicz et al. Phys. Rev. ST Accel. Beams 
7,044402 (2004). 

[9] K. Saeki et al. Phys. Rev. Lett. 42,8 (1979). 
[10] W. Fox et al. Phys. Rev. Lett. 101,255003 (2008). 
[11] G. Sam et al. Eur. Phys. J. D published online IOD: 

10. 1 140/epjd/e2009-001 15-8. 
[12] L. Romagnani et al. Laser Part. Beams 26,241 (2008). 



[13] M. J. Norman et al. Appl. Opt. 41,18 (2002). 

[14] J. D. Lindl et al. Phys. Plasmas 11,339 (2004). 

[15] R. A. Snavely et al. Phys. Rev. Lett. 85, 2945 (2000). 

[16] J. F. Dempsey et al. Med. Phys. 27,10 (2000). 

[17] M. Borghesi et al. Phys. Rev. Lett. 92,055003 (2004). 

[18] F. Califano, M. Lontano Phys. Rev. Lett. 95,245002 

(2005). 

[19] http://www.casinc.com/codes.html 
[20] B. Eliasson, P. K. Shukla Phys. Rev. Lett. 93,4 (2004). 
[21] M. V. Goldman et al. Phys. Rev. Lett. 99,145002 (2007). 
[22] D. Summers, R. M. Thorne J Geophys. Res. 97,16827 

(1992). 

[23] M. Borghesi et al. Appl. Phys. Lett. 82,1529 (2003). 
[24] K. Quinn et al. Phys. Rev. Lett. 102,194801 (2009). 



P. Guio et al. Nonlinear Proc. Geoph. 10,75 (2003). 

L. Romagnani et al. Phys. Rev. Lett. 101,025004 (2008). 



